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Abstract. In this "Invisible Universe" proceedings, we introduce the Dark Energy Universe Simulation Series (DEUSS) 
which aim at investigating the imprints of realistic dark energy models on cosmic structure formation. It represents the largest 
dynamical dark energy simulation suite to date in term of spatial dynamics. We first present the 3 realistic dark energy 
models (calibrated on latest SNIa and CMB data): ACDM, quintessence with Ratra-Peebles potential, and quintessence with 
Sugra potential. We then isolate various contributions for non-linear matter power spectra from a series of pre-DEUSS high- 
resolution simulations (130 million particles). Finally, we introduce DEUSS which consist in 9 Grand Challenge runs with 1 
billion particles each thus probing scales from 4 Gpc down to 3 kpc at z=0. Our goal is to make these simulations available 
to the community through the "Dark Energy Universe Virtual Observatory" (DEUVO), and the "Dark Energy Universe 
Simulations" (DEUS) consortium. 
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INTRODUCTION 

According to current Supernova la Hubble diagrams (Kowalski et al. [11]), and Cosmic Microwave Background 
anisotropies power spectrum (Komatsu et al. [10]), the universe energy budget is dominated by dark energy (about 75 
percents). However, the nature of dark energy is still unclear. What are the properties and origin of dark energy? How 
do the density and equation of state evolve across cosmic times? What is the spatial distribution? Several physical 
interpretations for the observed recent acceleration of the expansion of the universe have been proposed such as 
cosmological constant, backreaction (Buchert [4]), quintessence (Ratra & Peebles [18]) and deviations from general 
relativity (Alimi & Fuzfa [2]) among many others. In order to discriminate between these models one has to refine 
current observations at the homogeneous or linear level, or to find new observables. In this article, we focus on the 
latter idea by investigating the imprints of dark energy on cosmic structure formation using very-high resolution N- 
body simulations. We note that this article is complementary to Alimi et al. [1], Courtin et al. [7, 8], Rasera et al. 
[20]. 



REALISTIC QUINTESSENCE MODELS 

In this article, we consider three cosmologies with different properties: a quintessence model with a Ratra-Peebles 
potential (Ratra & Peebles [18]), a quintessence model with supergravity corrections (Brax & Martin [3]) and the 
standard ACDM model. We however expect our physical conclusions to apply on a much wider class of models. Most 
of the previous N-body simulations with dark energy assume ACDM cosmological parameters and even ACDM linear 
power spectrum for setting initial conditions (see however recent papers from Casarini et al. [5], Jennings et al. [9]). 
To overcome these issues, we create a set of two "realistic quintessence models" plus a reference ACDM WMAP-5yrs 
cosmology. We modify CAMB to take into account the fluctuations of quintessence in order to compute the linear 
matter power spectrum. We then perform a likelihood analysis of the UNION SN la Hubble diagram and WMAP- 




FIGURE 1. 68% and 95% confidence regions in the Q. m h 2 — a plane from the combined analysis of the UNION SN la Hubble 
diagram and WMAP-5yrs CMB data for the Ratra-Peebles quintessence models (left) and Sugra quintessence models (right). A 
cross mark (X) indicates our choice for realistic RP and Sugra model parameters, while some other RP and Sugra model parameters 
assumed in the literature (for simulations) are marked with other symbols. These figures come from Alimi et al. [1]. 

5yrs CMB data to select viable quintessence models as illustrated Fig. 1. We choose £l m ~ 0.23, a = 0.5, Ofe = 0.66 
for Ratra-Peebles quintessence models, D. m — 0.25, a = 1, o& = 0.73 for Sugra, £l m = 0.26, o& = 0.79 for ACDM. 
It is worth noting that quintessence models favor lower £l m and Og. Our three models are degenerated regarding to 
CMB and SNIa constraints and one needs additionnal constraints (such as from the non-linear regime) in order to 
discriminate them. 



HIGH-RESOLUTION N-BODY SIMULATIONS 

We run a series of 9 high resolutions simulations with 512 3 particles each, for the three "realistic" cosmologies de- 
scribed above and for three box lengthes (162, 648 and 1296 h _1 Mpc). Initial conditions are set using the MPGRAFIC 
software (Prunet et al. [17]). Particles evolution is computed with the Adptative Mesh Refinement (AMR) code RAM- 
SES (Teyssier [22], Rasera & Teyssier [19]). Finally, power spectra are computed using POWMES (Colombi et al. 
[6]). Resulting matter power spectra are shown in Fig. 2. From the power spectra plot (top-left), it can be seen that 
quintessence universes are less structured than ACDM ones at all redshifts. A large part of the discrepancies is im- 
puted to the different linear power spectra (lower Og but also quintessence-specific shape). In order to remove this 
linear contribution, we show (on the top-right plot) ratios of quintessence models power spectra (normalized by the 
linear one) to ACDM power spectrum (also normalized by the linear one) at z=0. As expected for small wavenum- 
bers (below 0.1 h/Mpc), the three power spectra superimpose. However, differences arise above 0.1 h/Mpc, which are 
mostly caused by non-linear amplification of the linear growth rate. This effect is well accounted for by the Smith 
et al. [21] fitting function. The bottom-left plot confirms this point for wavenumbers below 1 h/Mpc. However, some 
important discrepancies remain for k larger than 1 h/Mpc. Using a constant equation of state (0 = (Oq and the fit of 
McDonald et al. [15] doesn't solve the issue. It indicates that each dark energy model leaves a specific imprints on the 
high-k tail of the power spectrum, and that using constant equation of state approximation (Ma et al. [13], McDonald 
et al. [15], Lawrence et al. [12]) results in several percents errors on the power spectrum. The correlation with the 
linear growth history (bottom-right plot, see also Ma [14]) indicates that the non-linear power spectrum records the 
history of structure formation (Alimi et al. [1]). This has yet to be properly accounted for in analytical predictions 
since most of the current fits (Peacock & Dodds [16], Smith et al. [21]) only depend on the instantaneous value of the 
linear growth rate. 

GRAND CHALLENGE RUNS: DEUSS 

In order to gain further physical insight into dark energy's imprints on cosmic structure formation, one needs to 
investigate its signature at all scales: from galaxies inner profiles to Hubble volume cluster counts. To this aim, 
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FIGURE 2. In all the plots ACDM is the continuous line, Ratra-Peebles the dashed line, and Sugra the dot-dashed line. Top left 
panel: power spectrum at z=0 (upper triplet), 1 and 2.33 (lower tripplet) and associated linear power spectra (dotted lines). Top right 
panel: Ratio of power spectrum (normalized by the linear power spectrum) to the ACDM power spectrum. Bottom left panel: Ratio 
of the non-linear power spectrum (normalized by the Smith et al. [21] fit) relative to the ACDM one. The thin continuous line are 
the constant eos expectation for RPCDM (top) and SUCDM (bottom). Bottom right panel: Linear growth factor histories relative 
to the ACDM WMAP5 one (for comparison ACDM WMAPI and ACDM WMAPIII are shown as triple dot dashed line and long 
dashed line respectively). The deviations at high-k of the power spectra are correlated with the linear growth histories.These figures 
come from Alimi et al. [1] 



we perform the unprecedented Dark Energy Universe Simulation Series (DEUSS). It consists in 9 Grand Challenge 
simulations with 1 billion particles each (and up to 7 billions cells), resolving scales from 4 Gpc down to 3 kpc for 
the three considered cosmologies: ACDM, Sugra and Ratra-Peebles. Taking advantage of its Peano-Hilbert domain 
decomposition, the MPI code RAMSES has been run on 4096 Blue Gene/P cores at the IDRIS supercomputing center 
(last trimester 2008). It represents 5 millions mono-cpu hours (that is to say about 600 years) of computationnal time 
that have been obtained by the Horizon Project at LUTH. The run parameters are described in Table 1 : simulations 
box lengthes are 162 Mpc/h, 648 Mpc/h and 2592 Mpc/h. 

We backup three kinds of data during the runs. First, we save the masses, positions, velocities and identities of all 



TABLE 1. Simulation parameters for the DEUSS Grand Chal- 
lenge runs at the IDRIS supercomputing center. Each simulation 
contains 1024 3 particles with 1024 3 coarse-grid cells, 6 refinement 
levels and has been evolved down to z = 0. The various entries report 
the values of the simulation initial redshift Zi, particle mass m p and 
spatial resolution A x . 



Parameters 


162h" 1 Mpc 


648 fr'Mpc 


2592h~ 1 Mpc 




ACDM W5 

Zi 

m p (h- x M ) 
A^h^kpc) 


137 
2.86 x 10 8 
2.47 


93 

1.83 x 10 10 
9.89 


56 

1.17 x 10 12 
39.6 


Ratra-Peebles 

Zi 

m p (\T x M©) 
A^h-ikpc) 


118 
2.53 x 10 8 
2.47 


81 

1.62 x 10 10 
9.89 


50 

1.04 x 10 12 
39.6 


Sugra 

Zi 

m p (hr x M Q ) 
A^h-ikpc) 


134 
2.75 x 10 8 
2.47 


92 

1.76 x 10 10 
9.89 


56 

1.13 x 10 12 
39.6 



particles for 24 snapshots for each simulation (expansion factor a = 0.025, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 
0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.729, 0.75, 0.8, 0.823, 0.85, 0.9, 0.95 and 1), that is to say a total of 216 snapshots. 
Second, we construct 9 lightcones around an observer at the center of the simulation by dumping particles from the 
light-travel distant slice at each coarse time step. For the largest box, the lightcones cover An steradian up to z=l 
whereas for the smallest volume, the opening angle is 10 degrees by 10 degrees. Third, we save 6 samples of particles 
at every coarse time step in a parallelepiped of size 1/16 x 1/16 x 1 times the full volume of the simulations. The 
largest number of coarse time steps is 1419. The total amount of data is about 40 To. 

Unfortunately, the output files are organized along the Peano-Hilbert curve which could be non-trivial to manipulate. 
They also contain a lot of double-precision and not-useful information (except for restarting the run). We therefore run 
a parallel Sheer in order to split the data in 8 3 cubic regions or "fields" and to keep only useful simple precision data. 
We also develop and run a parallel Friend-of-Friend halo finder (linking length b=0.2) on 512 processes. We detect 
all halos above 100 particles in half of the outputs (the analysis of the others is on-going). We then sort the particles 
halo per halo and backup these 512 "halo" files. Overall halo masses range from 3 x 10 10 Mpc/h to 8 x 10 15 Mpc/h, 
there are up to 500 000 halos per snapshot and up to 3 millions particles per halos. The total number of analyzed halos 
is about 20 millions (20 other millions will be detected soon). With this two kinds of files (halos and fields), most of 
subsequent post-processing can be serial. Moreover, data are user-friendly to manipulate. 

As an illustration of the results we show a composite view of the density field for the 3 simulations of box length 
162 Mpc/h at z=0 with color coding for each cosmology: ACDM (red), Sugra (green) and Ratra-Peebles (blue). On 
large scale, the map is essentially black and white indicating little difference between cosmologies whereas on smaller 
scales (40 and 10 Mpc/h) the differences between the models show up (multi-color filaments and halos). Discrepancies 
are therefore amplified by the non-linear regime of structure formation and can help to break degeneracies between 
cosmologies. 

Our plan is to share these data with the "dark energy" community. To this aim we are starting the DEUS (Dark 
Energy Universe Simulations) consortium. The goal is to gather observers and theoreticians in order to analyse the 
data, prepare future dark energy observational missions (weak lensing, strong lensing, BAO, cluster counts, etc.), 
improve analytical predictions (power spectrum, mass functions, PDF, density profiles, etc.), and better understand the 
imprints of dark energy on structure formation. We are also building a database to diffuse the results. The database 
should contain all the properties of all halos and fields from the different outputs. The strength of this Dark Energy 
Universe Virtual Observatory (DEUVO) is that it follows the IVOA data model to allow future interoperability. A 
website is also currently under development (http://www.deus-consortium.org). 




FIGURE 3. Projected density maps from the Dark Energy Universe Simulation Series (DEUSS) at z=0 for three cosmologies: 
ACDM (red), Sugra (green) and Ratra-Peebles (blue). The spatial dynamics is illustrated by this zooming sequence from the full box 
scale (162 Mpc/h) onto one halo (at 10 Mpc/h scale) among the 300000 well-resolved halos in the simulation volume. At large scale, 
the picture is essentially black and white because differences between cosmologies are small. However, the imprints of dark energy 
(multicolor pictures) are clearly visible on non-linear scales. Overall, DEUSS consist in 9 Grand Challenge simulations with 1 
billion particles resolving scale from 4 Gpc to 3 kpc. It has been run on 4096 Blue Gene/P cores at the IDRIS supercomputing center, 
using 5000000 mono-cpu hours thanks to Horizon Project at LUTH. This figure comes from http://www.idris.fr/docs/docu/projets- 
Babel/DEUSS/CR-projet-DEUSS.html. 

CONCLUSION 

This article can be summarized as follow. 

1. Unlike most of the previous numerical simulations, we consider three realistic dark energy models. These models 
are calibrated on latest SNIa and CMB data and take into account quintessence clustering in the linear regime. As 
a result, £2„, and <7g are lower for quintessence cosmologies. 

2. We run a series of pre-DEUSS high-resolution runs (130 millions particles). We distinguish three dark energy 
signatures on the non-linear power-spectra: linear normalisation and shape, non-linear amplification, and a record 
of the past history of cosmic structure formation (which is not accounted for in most of current fits). 

3. We run the unprecedented DEUSS Grand Challenge runs in order to probe structure formation at higher resolution 
and on larger scale than before (ie. from 3 kpc to 4 Gpc) for the three realistic dark energy models. The simulations 
evolve 1 billion particles and up to 7 billions AMR cells. Our data consist in 216 snapshots, 9 lightcones and 6 
samples. We have already analyzed half of the snapshots and produced a catalog of 20 millions halos and 60 
thousands cubic fields. These user-friendly data represents a unique way for observers and theoreticians from the 



(newly-created) DEUS consortium to seek for new observables using the DEUVO database (in construction). It 
will help us to break degeneracies between dark energy models and to improve our understanding of the role of 
dark energy on cosmic structure formation. 
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